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A very simple way to obtain comparable baryon and DM densities in the early Universe is through 
their contemporary production from the out-of-equilibrium decay of a mother particle, if both pop¬ 
ulations are suppressed by comparably small numbers, i.e. the CP violation in the decay and the 
branching fraction respectively. We present a detailed study of this kind of scenario in the context 
of a R-parity violating realization of the MSSM in which the baryon asymmetry and the gravitino 
Dark Matter are produced by the decay of a Bino. The implementation of this simple picture in a 
realistic particle framework results, however, quite involving, due to the non trivial determination 
of the abundance of the decaying Bino, as well as due to the impact of wash-out processes and of 
additional sources both for the baryon asymmetry and the DM relic density. In order to achieve a 
quantitative determination of the baryon and Dark Matter abundances, we have implemented and 
solved a system of coupled Boltzmann equations for the particle species involved in their genera¬ 
tion, including all the relevant processes. In the most simple, but still general, limit, in which the 
processes determining the abundance and the decay rate of the Bino are mediated by degenerate 
right-handed squarks, the correct values of the DM and baryon relic densities are achieved for a 
Bino mass between 50 and 100 TeV, Gluino NLSP mass in the range 15-60 TeV and a gravitino 
mass between 100 GeV and few TeV. These high masses are unfortunately beyond the kinematical 
reach of LHC. On the contrary, an antiproton signal from the decays of the gravitino LSP might be 
within the sensibility of AMS-02 and gamma-ray telescopes. 
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I. INTRODUCTION 

The origin of the Dark Matter (DM) component of the Universe and of the baryon asymmetry 
are two compelling puzzles of modern particle physics and cosmology. Conventionally, different 
and unrelated mechanisms are considered for the generation of these two quantities. Indeed, the 
generation of the correct baryon asymmetry is a rather difficult task to achieve and requires definite 
conditions [T] to occur, mainly consisting in efficient B and CP violating processes occurring outside 
from thermal equilibrium. On the contrary the correct DM relic density can be generated by a very 
broad variety of mechanisms, also compatible with thermal relic DM particles, like the popular 
WIMP paradigm. 

A common generation mechanism for the DM and the baryon asymmetry is nonetheless a very 
intriguing possibility, also motivated by the similarity of the values of the two relic densities, being 
indeed Uab/Hdm ~ 0.2. 

The most simple way to connect the baryon and the DM abundances is to assume also for the DM a 
generation through asymmetry. In the simplest realization of asymmetric DM models (see e.g. [2]) 
the ratio Uab /^dm simply corresponds to the ratio between the mass of the proton and the mass 
of the DM. 

A viable alternative is however represented by the possibility of linking the generation of the baryon 
asymmetry to the popular WIMP mechanism. In this kind of scenarios the baryon asymmetry is 
produced after the chemical decoupling of a thermal relic, through B and CP violating annihila¬ 
tions PH or decays mm, or even of the DM itself QUIiP- 

In a similar spirit, a rather simple mechanism allowing to achieve the contemporary production of 
the DM and of the baryon asymmetry has been proposed in m■ Here these two quantities are 
contemporary generated by the out-of-equilibrium decay of a mother particle. The ratio between 
the DM and baryon density is expressed in terms of two analogously suppressed quantities, namely 
the CP-asymmetry and the branching ratio of decay of the mother particle into DM, and can be 
accommodated to be of the correct value, irrespectively of the initial abundance of the decaying 
particle, through a suitable choice of the parameters of the underlying particle theory. 

Although simple and elegant, this idea may be rather difficult to be implemented in concrete particle 
physics frameworks. For example, although the ratio between the baryon and DM abundances is 
independent from the one of the mother particle, this is not the case for the individual expectations 
of these two quantities. In order to match their rather precise experimental determinations m, 
a similarly precise determination of the abundance of the mother particle is required. This is 
in general a not trivial task since it is determined by many different processes. In particular 
additional states of the underlying particle theory might play a relevant role through coannihilation 
effects. Furthermore, in presence of extra new particles, with respect to the mother particle and the 
DM, additional sources of baryon asymmetry and DM can be present, spoiling the simple picture 
discussed above. Analogously crucial is finally the determination of the impact of possible wash-out 
processes, namely the processes capable of depleting a possibly generated baryon asymmetry. 

In order to properly deal with these issues, a detailed numerical treatment, relying on suitable 
Boltzmann equations, is mandatory. 

In this work we will investigate a definite case of study, being a R-parity violating realization of 
the MSSM with gravitino DM. SUSY models are a rather good playground for the scenario under 
consideration. Indeed, in absence of R-parity the SUSY superpotential features automatically 
sources of Baryon and Lepton number violation. At the same time, thanks to its Planck suppressed 
interactions, the gravitino DM remains stable on cosmological scales even in absence of symmetries 
forbidding its decay. In addition the decay branching fractions of supersymmetric particles into 
gravitino are as well Planck suppressed, thus not preventing an efficient generation of the baryon 
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asymmetry from the out-of-equilibrium decay of a supersymmetric state. 

The mother particle is instead a Bino-like neutralino which generates the baryon asymmetry and 
the DM by late-time decays occurring after its chemical freeze-out. These two quantities are sub¬ 
stantially determined, with mild assumptions on the cosmological history, by the underlying particle 
physics framework, in particular by the structure of the Supersymmetric spectrum. The require¬ 
ment of the correct baryon and DM abundances can be translated into predictions on the particle 
content of the theory. These predictions can be tested by collider experiments if the particles 
involved in their generation are within their production reach. 

We have determined the baryon and DM relic densities, Dab and Ddm through a system of coupled 
Boltzmann equations tracing the time evolutions of all the particle species involved in the generation 
of these two quantities. These are the Bino and the other two gauginos, the Gluino and the 
Wino, while scalar superpartners and the Higgsinos should be set, as clarified below, to very high 
scales such that they do not directly enter the system of equations as particle species, but only as 
mediators of the interactions of the gauginos. The system finally includes two additional equations, 
respectively for the baryon and the DM abundances. This kind of system has been solved as function 
of the relevant supersymmetric parameters. In order to provide a more clear understanding we will 
complement, where possible, our numerical treatment with analytical expressions for the relevant 
quantities. We will envisage, in particular, the dependence of the rates of the relevant processes, as 
well as the CP-asymmetry, on the flavor structure of the theory. 

The paper is organized as follows. After a brief review of the general idea of the contemporary 
generation of the DM and baryon densities from decay of a thermal relic, we will present in section 
III its MSSM realization. We will then present in section IV some analytical estimates of the relevant 
quantities. Section V will be instead dedicated to the numerical treatment and the quantitative 
determination of the parameters space compatible with the experimental expectation of the baryon 
asymmetry and of the DM relic density. Before stating our conclusions, we will finally briefly 
mention in section V the possible detection prospects relative to our setup. 


II. DARK MATTER AND BARYON PRODUCTION FROM OUT-OF-EQUILIBRIUM 

DECAY 

A simple and elegant way to achieve the contemporary production of the baryon asymmetry and 
of the Dark Matter is, as proposed in |14j . by out-of-equilibrium decay of a state X, featuring B 
and CP violating interactions, and thus capable, according the Sakharov conditions, of generating 
a baryon asymmetry. The resulting baryon density can be schematically expressed as: 

TTL _ 

Dab = Cabecp —— BR(X —► b , b) D.y (1) 

m x 

where m p is the mass of the proton, ecp is the CP asymmetry: 

_ r(x ->■ b) - r(x -> 1) 

CP T(A -a b) + T(A -a b) ' U ' 

and Da is the initial abundance of the state A'. The factor £ab encodes the effects of the sphaleron 
processes, as well possible wash-out and entropy dilution effects. The field X features as well an 
additional decay channel, not necessarily B-violating, into DM such that its relic density is given 
by an analogous expression as above: 

Dbm = BR (A -A DM + anything) D x (3) 


m x 
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Expressions (JTJ) and ^ both feature suppression factors. For example the baryon density is sup¬ 
pressed by the ratio m p /mx- In addition, in most realistic particle frameworks, the CP asymmetry 
ecp is a suppressed quantity. At the same time it is reasonable to expect that the branching ratio 
of decay of the X particle into DM is as well suppressed in order to do not dangerously affect 
the baryon production. As a consequence in order to account for the experimental expectations 
of Ddm and Dab a rather high value of the initial D\- is needed. In addition, under the assump¬ 
tion, performed in this work, that the initial abundance of X is generated similarly to the WIMP 
mechanism, we need to require for it a sufficiently long lifetime such that it decays after chemical 
freeze-out. 

Interestingly, the ratio of the two densities is independent from the one of the X state, being: 

Dab ^ m p BR(X — » b , b ) £ab , 

D dm 6CP rn DM BR ( X -> DM + anything) ’ £dm 

and its expected value ~ 0.2 is achieved by a suitable choice of the DM mass and of the param¬ 
eters determining the CP asymmetry and the two branching ratios. In this work we will embed 
this mechanism in a supersymmetric framework with gravitino DM while the decaying state X is 
represented by a Bino-like neutralino. 


III. MSSM REALIZATION 

In our investigation of the possibility of contemporary production of the baryon asymmetry and 
of gravitino DM we will focus on a MSSM realization with R-parity broken only by the operator 
A U C D C D C which provides the breaking of the baryon number avoiding at the same time strong 
constraints from the stability of the proton Q being lepton number violating operators absent. 
The mother particle is chosen to be a Bino. As will be shown below, the rates of the processes 
governing its abundance and lifetime are set by the mass scales of the scalar superpartners and of 
the Higgsinos [^] The requirement of overabundance and long lifetime of the Bino can be met for 
Supersymmetric spectra like the ones proposed in nuns featuring a strong mass hierarchy between 
the three gauginos, namely the Bino, the Gluino and the Wino, and the scalar superpartners as well 
as the Higgsinos. The Bino is not, however, the Next-to-Light Supersymmetric particle (NLSP). The 
CP-asymmetry is created from the interference between tree-level and one loop processes involving 
a quark, a squark and another gaugino [5j. According the Nanopolous-Weinberg theorem njang- 
1 2D, this asymmetry is not null (at this order in perturbation theory) only if at least one between 
the squark and the gaugino running in the loop is lighter than the Bino. The case of a squark 
NLSP is however not feasible in a MSSM setup since a light-squark would enhance the annihilation 
and decay rate of the Bino making it not enough abundant and long-lived. A viable scenario can 
be obtained, for example, by extending the MSSM with an additional singlet playing the role of 
mother particle mm- 

In the scenario considered the lightest gaugino is the Gluino. It is possible to realize, alternatively, 
a leptogenesis scenario by considering the operator A QLD C , rather than the R-violating one, and 
considering the Wino as lightest particle, apart the gravitino DM. The Gluino and the Wino are 


1 The proton can actually decay into a gravitino, if kinematically possible, even in presence of only the A cou¬ 
pling m ■ As will be seen at the end of the paper the favored region of the parameter space will feature a gravitino 
much heavier than the proton, such that this kind of decay is forbidden. 

2 We can anticipate the most of the relevant processes occur before the Electroweak Phase Transition (EW) tem¬ 
perature. As a consequence the spectrum of the electroweakly interacting fermionic superpartners consists of two 
Majorana fermions, the Bino and the Wino, and two Dirac fermions, the Higgsinos. 
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instead not good candidates for the generation of the baryon asymmetry since they feature very 
efficient annihilation processes into gauge bosons, with rates depending on their same masses, which 
make their abundances too suppressed to generate sizable amounts of baryons and DM. 

The processes responsible for the generation of the baryon asymmetry can be described by the 
following effective lagrangian [22], where, in agreement with the discussion above, the scalars and 
the Higgsinos have been integrated out (for simplicity we are omitting the mass terms and, from 
now on, indicate by A, rather than, as conventional, A , the RPV coupling.): 


C = —2e' 


aft 7 


+A kij B 

+Afcij | B 


{*!« [B (< 
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(5) 


where: 
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where T^ Lai are 6x3 matrices defined by: 


Qa — r LaiQLi T C RctiQRi 


( 6 ) 

(7) 


where q a , a = 1, ■ • • ,6 are the squark mass eigenstates while qL,m , i = 1, ■ ■ • ,3 are left-handed 
and right-handed squarks. The quantities g are instead defined as: 

gf = -V2 9l ( Q f - T 3 ) e** ff| R = V2 9l Q f e^R 

g^ = -V2g 3 e^o 5 | R = (8) 

where (j)^ and (j )^ are CP -violating phases (see below for more details). 

The lagrangian pi) gives rise to CP violating decays of the gauginos into three SM fermions, e.g. 
B —> udd, as well as CP violating 2 —> 2 scatterings. The baryon asymmetry arises in the decay 
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(and annihilations) of the Bino, through the interference of tree-level diagrams, generated by the 
interactions proportional to the B-violating couplings A, with loop level diagrams obtained by 
combining these operators with the effective B — G interactions reported in the last lines of (15j) . 
Inverse decays and 2 —> 2 scatterings involving all the gauginos, in particular the Gluino, represent 
the main wash-out processes which guarantee that no baryon asymmetry is created in thermal 
equilibrium. Moreover, B-violating single annihilation processes Bd(u ) —> d(u)d play a prominent 
role in determining the abundance of the Bino, together with the Bino-Gluino coannihilations. 
The very last line of |5| is not associated to any process involved in the generation of the baryon 
asymmetry but triggers the pair annihilation BB —> HH* which is also relevant for determining 
the abundance of the Bino and, consequently, the one of DM and baryons. As evident the effective 
coupling depends on the combination sin /3 cos /3. In the analytical expressions provided below we 
will implicitly assume the limit tan/3 —> 1, in order to guarantee the correct EWSB and avoid 
tensions with the determination of the Higgs mass mmmi given the high scalar mass scale (See 
however the recent analysis [24]). 

The expressions of the relevant interactions rates are in general very complicated, in particular 
because of a non-trivial interplay of the flavor structure, which is substantially free due to the very 
high scale of the scalar masses. However, as clarified in the next section, we can present our results, 
without loss of generality, in a simplified limit in which effects of flavor violation are neglected and 
the relevant interactions are mediated by only down-type right-handed squarks. This choice allows 
for simpler computations, since the number of possible processes is reduced; at the same time it 
guarantees the presence of all the possible topologies of diagrams responsible of the generation of 
the baryon asymmetry. For analogous reasons we have not reported the interactions of the sleptons 
which are assumed to be decoupled. In general one should consider analogous operators as the 
ones reported in eq. <© also for the Wino. These however, as further discussed later on, do not 
significantly contribute to the generation of the baryon asymmetry and thus have been omitted for 
simplicity. 

The DM candidate in this scenario is the gravitino. Although not exactly stable in an R-parity 
violating scenario, the Planck suppression of its interactions guarantees a lifetime largely exceeding 
the one of the Universe, even for 0(1) values of the RPV couplings [25] [26]. The Bino (as well 
as the other superpartners) features a decay channel into the gravitino LSP with a Planck scale 
suppressed branching fraction. The generation of the baryon asymmetry is then insensitive to this 
decay channel which can nonetheless produce a sizable amount of Dark Matter since the suppressed 
branching ratio can be compensated by the overabundance of the decaying mother particle. The 
out-of-equilibrium decay of the other gauginos does not, instead, efficiently produce DM in view of 
their suppressed relic abundance. The Wino and the Gluino, as well as the Bino itself, can anyway 
copiously produce DM, by freeze-in, at early epochs while they are still in thermal equilibrium. 

In order to have a proper description of the generation of the baryon asymmetry, as well as the 
DM relic density one must rely on a solution of system of coupled Boltzmann equations. Before 
illustrating it, we will anyway provide some analytical approximations in order to provide a better 
understanding of the results, in particular the implications for the Supersymmetric spectrum. 


IV. ANALYTICAL RESULTS 

We will present in the following analytical expressions for both the baryon and the DM relic den¬ 
sities. These expressions are strictly valid only in definite regions of the parameter space while a 
systematic investigation requires the solution of a system of coupled Boltzmann equations, like the 
one presented in the next section. For greater clearness we will discuss separately, in the next two 
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subsections, the generation of the baryon and of the DM abundances. 


A. Generation of baryon asymmetry 

As already mentioned, there are actually two sources for the baryon asymmetry, namely the B 
and CP violating decays of the Bino as well as 2 —> 2 scattering processes. In the first case the 
asymmetry is originated by the three-body decays, B —> udd ( udd ). In the case that these processes 
are mediated by only right-handed d-squarks (see discussion below) the relevant tree-level and one- 
loop diagrams are shown, respectively, in fig. 0 and ([2]) (see also pH HI]). 




(a) Diagram T1 


(b) Diagram T2 


FIG. 1: Diagrams contributing, at the tree level, to the B violating decay of the Bino in the case of mediation 
from only d-squarks. 


A CP asymmetry is generated as well by annihilation processes like, e.g., Bd —> ud + CP conjugate. 
The relevant diagrams are obtained, by crossing symmetry, from the ones shown in figs. 0-0- 
Assuming distinct timescales for the B -violating scatterings and decays, the total CP asymmetry 
ecp can be expressed as: 


£CP 


Ar dec 

I tot .dec 


Ar a 


(9) 


where Ar ann and Ar dec are, respectively, the differences between the rates B violating scattering 
and decay processes and their CP conjugates. The corresponding expressions are: 
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(b) Diagram L2 




(c) Diagram L3 


(d) Diagram L4 


FIG. 2: Diagrams contributing, at the loop level, to the B violating decay of the Bino in the case of mediation 
from only d-squarks. The CP-asymmetry is generated by the interference with tree-level diagrams reported in 


fig-Q- 
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(14) 


while Ai 2 , Bi 2 ,C\ 2 , Di 2 are numerical coefficients which can be determined from the expressions 
provided in the appendices. 

A non-null CP asymmetry originates from a non trivial combination of the phases coming from from 
the (Majorana) gaugino masses, encoded in the couplings g^ and gg, the squark mixing matrix and, 
possibly, the RPV couplings A. We notice in particular that the terms proportional to the coefficients 
A\ 2 and Ci ,2 are different from zero only in presence flavor violation since the combinations between 
the gauge couplings are automatically real and the phases in the RPV couplings would as well cancel 
in this limit. In absence of flavor violation the CP violation arises from the differences of phases 
contained in the Maiorana masses of the Bino and the Gluino, which behave as effective B- violating 
terms 1312110 In this case however the CP asymmetry is suppressed by the ratio as well 

as by the kinematical factor f 2 . As already mentioned we are considering a regime in which only 
right-handed d-type squarks contribute to the processes of interest. As clarified in the appendix 
additional contributions are originated by similar diagrams in which up-type right-handed squarks 
are exchanged. However the eventual increase of the CP asymmetry does not necessarily imply an 
increase of the baryon abundance. Indeed there is a tight relation between the processes governing 
the generation of CP asymmetry with the ones governing the abundance of the Bino as well as 
wash-out processes. In general an increase of the CP asymmetry is connected with an enhancement 
of the depletion rates of the Bino and of the baryon asymmetry itself and one has then to find a 
balance between the two effects. This provides a further indication that any analytical treatment 
should be complemented by the numerical solution of suitable Boltzmann equations. 

The suppression rriQ / m^ can be also avoided in presence of mixing between left and right-handed 
squarks, which would make to arise analogous terms as the first in eq. PPI ). We remind however 
that the size of left-right mixing depends on the ratio Xf /m§ where: 


X f = m f ( A f - nqp) 

tan /3 for d-type squarks 
cot P for u-type squarks 


qy = 


(15) 


where m/ is the mass of the SM fermionic partner of the squark. This mixing is thus heavily 
suppressed as the ratio rrif/mq a with the only possible exception of the top squark, where the 
ratio m t /mq a might be balanced by taking 1 (The Af terms can differ at most by a 

0(1) factor from rriq in order to avoid breaking of the color.). As clarified below, an efficient 
production of the baryon asymmetry requires rriq a > 10® GeV. For such values we can achieve 
values of X t /m^ a ~ lO - ^ 1 - 1 which do not induce sensitive variations of the total CP asymmetry 
with respect to the simplified regime we are considering. 

Even in presence of flavor violation the contribution from the coefficients A\ 2 and C\ 2 is limited 
since the combination between the flavor matrices is GIM suppressed and its imaginary part is zero 


3 We are encoding the CP-phases in the vertices gaugino-quark-squark while the Majorana masses of the gauginos 
are assumed to be real. This configuration can be obtained through a suitable rotation of the superfields. 
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in the limit of degenerate squarks. Even for non-degenerate squarks it is possible to achieve at most 
0(1) variations of the CP-asymmetry with respect to the flavor universal scenario. The remark 
again that the impact of this variation is not trivial to identify at the analytical level because of 
the non trivial interplay with the wash-out processes and the ones responsible of the abundance of 
the Bino. We will thus postpone further discussion of this point to the section dedicated to the 
numerical analysis. 

We finally notice that the functions fi and / 2 in (10) as well as /AEj and IA£ 2 in (121 make the 
asymmetry zero, consistently with Nanopolous-Weinberg theorem, if t/Iq > mg. 

In agreement with the discussion above, without loss of generality, we will present our results in 
the limit of absence of flavor violation and degenerate squark masses m,q a = mo- In this limit the 
expressions above simplify to: 


AT dec = 


AT — 

—* 1 arm 


m%,mr 



Im [e 2l?i ] 


mX m 


, —— ) AngTm [e 2 


m 


= ~f ( 16 ) 


where (f> = <f>Q — </g. As further simplification we will assume that all the couplings A kij (apart the 
ones set to zero by the asymmetry of the ij indices) are equal to the same value A. 

A sizable asymmetry from 2 —» 2 scatterings can be created only for freeze-out temperatures of the 
Bino very close to its mass. For lower temperatures, indeed, it results drastically reduced by the 
Boltzmann suppression in A?rg |21| . On the other hand, at high temperature, wash-out processes 
are still active and tend again to reduce the contribution to the asymmetry. As shown in j5] the 
correct amount of baryon asymmetry from 2 —» 2 scatterings can arise only from a very restricted 
range of values of the relevant parameters. On the contrary out-of-equilibrium decay can lead to 
a very efficient baryon production since it occurs at later time and, as clarified in the following, 
can evade wash-out effects if the Bino is enough long-lived. As also confirmed by our numerical 
investigation the decay of the Bino accounts for substantially the total amount of baryon density 
in all the viable regions of the parameter space. The baryon density reduces to 0 : 

b = £a b —-ecp fig 4 ' 00 (17) 

TOg 


The parameter £a b can be decomposed as the product £a b = £sp£w.o.£s- £ sp represents the effects of 
the sphaleron processes and can be set to 28/79 or 1 depending on whether the Bino decays before or 
after the temperature of electroweak phase transition, set to Tew = 140 GeV. £ w .o. and £ s represent 
instead the possible reduction of the baryon abundance due to wash-out effects while £ s is related 
to possible entropy dilution effects. An analytical estimate of the latter is provided at the end of 
this subsection. We have instead no analytical estimation for £ wo .. We can nonetheless identify, as 
explained below, two limit regimes, namely the case £ w .o. •C 1, corresponding to negligible baryon 
abundance, and £ w 0 = 1, for which a viable phenomenology is instead achievable. 

In agreement with what stated above, the CP asymmetry ecp is given by: 


AT dec 
ecp = -p- 

J- tot 


(18) 


with: 


Ttot = r ^B —> udd + udd^j + (j3 —> Gdd'j + T (B ^ 3 / 2 + 


( 19 ) 
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where A' represents all the possible SM final states accompanying the gravitino and: 


r (B —> udd + udd ) = 


A 2 ai 
167t 2 mg 




(B Gff) 


( 20 ) 

( 21 ) 


r(s->4 /2 + x) 


1 m l 

487 t mg^ 2 M|j 


( 22 ) 


are, respectively, the tree-level B-violating decay rate of the Bino in three SM fermions and of 
the B-conserving channel into the gluino and a pair of d-quarks, and finally the decay rate into 
any final state with gravitino, responsible of DM production. Mpi is the reduced Planck mass 
Mpi = 2.43 x 10 18 GeV. The last decay channel does not affect the baryogenesis mechanism in view 
of its very suppressed branching ratio: 


Br 


(B -> 4/2 + Aj w 5.7 x lO" 10 ( 1 + 


6A 2 


The CP asymmetry is then given by: 


£cp = ? Jm [e 2 


-a s (1+ 


TTOsN 

6A 2 / 


/ m 3/2 \ 

_2 / mo \ 

V lGeV / 

V 10 6 GeV J 



(23) 


(24) 


From now on we will take the value of the phase giving maximal ecp and assume that Im [e 2 *^] = 1. 
We notice that for A > WBiE"' the CP asymmetry is substantially independent from the amount of 
R-parity violation. For lower values it instead decreases as A 2 . The CP asymmetry is suppressed 
by the ratio as well as by the kinematic function / 2 . In order to achieve the correct 

baryon abundance this suppression should be compensated by a sufficiently high initial abundance 
of the Bino, which is set by its annihilation processes. 

These are described by thermally averaged cross-sections which can schematically be expressed asQ 


1 

W v ) (XiXj “► XiXk) = &rr 2 2 / m .x fmi \ / ds PijWijK i 

8Tmfm 2 A 2 (^) A 2 (^f) J^+rntf 


T 


w » = J dSl |Mp 


Pij = 


\]S - (rrii - iTi j )^ s - {rrii +rrijf 


2^s 


(25) 


4 The extrema of integration in principle exceeds the energy scales for which the effective description |S]) is valid. 
However, as will explained in the next subsection, in order to have a cosmologically viable scenario, we need 
to assume a low reheating temperature such Tr < mo- As a consequence, all the rates will be the computed at 
temperatures such that neglecting the momentum dependence of the propagators of the squarks and the Higgsinos, 
does not produce sensitive variations in the results. 
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The possible annihilation processes include, first of all, conventional pair annihilations; in our sce¬ 
nario the dominant ones are into two Higgses or into two SM fermion final states. The corresponding 
cross-sections are: 


{(TV 




A(x) =L dzz{z2 - fa2) Ki[z) 


(av) (bB -> 99 ) = 


m 


o 


Y K 3 (x)\ 2 

( Ki (x) A 2 " 

\k 2 (x)J 

\K2(X)J _ 


We have then coannihilation m processes with the two other gauginos: 


(crv) yBG —>■ uuj + (av) yBG —> daj = 
(av) (BW -> HH*) = 


_ 167 raia s m\ f^ K 4 (x) ^ 


27 


J K 2 {x) 


aiazn / TOg m^r 

32 y 2 \ T m.g 


(26) 


B(x,y) = 


x i y 2 K 2 (x)I\ 2 (yx 


- j dz (z 2 -±x 2 (l+y) 2 f /2 (z 2 -Ax 2 (l-y) 2 ) 1,2 K 1 (z) 


(27) 


We remark that a sizable contribution from coannihilations with at least one gaugino is unavoidable 
in our scenario since, in order to have a non-zero baryon asymmetry, the presence of a lighter 
gaugino with respect to the Bino is mandatory. Contrary to conventional WIMPs, we have to take 
into account also single annihilation processes, both RPC and RPV, with a SM fermion as second 
initial state. The relevant cross-sections are: 


) ^ Bu —>■ Gu'j + (av) (j3d — > Gd^j = 

) (^Buk —> didj'j + (av) (j3di —> Ukdjj 


47 ro:ia s rn-~ / K 4 (x) 
27 Wg \ K 2 (x) 


+ 1 


aqA 2 


nrif 


5^4(z) 

5 K 2 (x) +1 


(28) 


The relative contributions of the various annihilation channels, expressed in the form T ann /H where 
H is the Hubble expansion parameter and T ann = (av)n B y where X = B for pair annihilation pro¬ 
cesses, X = G,W for coannihilations, and X = q for single annihilations, are shown in fig. ^ . We 
have considered there four assignations of the set (nig, mo, y) while we have fixed the remaining 
parameters as itiq = 0.35 m^, m^, = 5wg, A = 0.2. The pair annihilation cross-section, in par¬ 
ticular the HH* channel, dominates for lower masses of the Bino and small hierarchy between nio 
and /i while, once increasing these quantities, single annihilation processes are the most important 
in determining the abundance and decoupling time of the Bino. Single annihilations are also in 
general dominant for low value of mo- We notice from the second and third panel of fig. Q, with 
mo set, respectively, to 10 6 and 10 6 ’ 5 GeV, that single annihilation processes determine a very late 
(even more than conventional WIMP scenarios) chemical decoupling of the Bino, for which a very 
suppressed relic abundance is expected. The results shown in fig. §} thus provide a first quali¬ 
tative indication that very high values of the scale mo are required to generate a sizable baryon 
abundance. A quantitative determination of the abundance of the Bino necessarily relies on the 
numerical solution of Boltzmann equations, illustrated in the next section, in particular because 
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FIG. 3: Annihilation rates, normalized with the Hubble expansion factor H, for the channels reported in the 
plot, of the Bino, for four assignations of reported in the panels and rriQ = 0.35 m ^, rn = 5 m ^, 


A = 0.2. 


single annihilation processes can induce deviations from the conventional WIMP scenarios. Indeed 
an analytical estimate of the baryon adundance is given by: 


Yg(X f) = M(x f) 


M(x j) 
y b( x 0 


1 “I 


!>}l Y q ,e 


(M(xi) -M(x f )) 


M(x) = exp -(av)\Y qt 


1 §M pi 


(29) 


where ( av) p and (<jv)i represent, respectively, the sum of the thermally averaged pair (including 
coannihilations |27| 1 and single annihilation cross-sections. Y q , eq = n g eq /T represents the yield of 
the quarks (constant in the relativistic limit). [^] Xi and X{ represent, respectively, an initial time, 


5 In writing eq. (|29|) we have neglected the time dependence of the annihilation cross-sections, in order to provide a 
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which can determined through an analogous procedure as presented in [28], and a final time which 
can set to be the decay time scale of the Bino, as defined below. In the limit a(crv)\/x <C 1 and for 
late enough decays (such that the first term in the parenthesis can be neglected) it is possible to 
recover the conventional WIMP behaviour, Y (x f) oc . The relic density, in this limit, is well 
approximated by the well known formula [28] : 


3.9 x lOVo. GeV" 1 
gl /2 Mpi (av(X{, 0 ,)) p 


(30) 


where the effective thermally averaged pair annihilation cross-section is computed at X[, Q = ^—- 
with Tf. 0 being the freeze-out temperature. In the limit in which the dominant annihilation channel 
is the one into HH* the Bino abundance is given by the rather simple expression: 


nr-^oo 

lL B 


4.1 x 10 9 


10 8 GeV 


2 Xj.o. 

A( x i.o.) 


(31) 


We can also expect that for high enough values of the scales mo and /z, the consequent suppression 
of the annihilation cross-section leads to a relativistic decoupling of the Bino. In such a case its 
relic abundance would be even larger [25]: 


or-^oo _ 7 o x in 10 9 b ( \ 

“ B ,rel -7.8x10 g ^ X{ ^[ lTeY ) 


(32) 


By using eq. (31) and (24), and setting £ sp = £ w .o. = £s — 1, we can write the baryon abundance 
as: 


n AB h 2 « 3.3 x 10 


-2 gf.o. ( m B \ ( f ( 

A(x{. 0 .) vITeV/ \rnfj) 2 ym| J 




10 3 / 2 to o / \7 ra, 


^ ( 1 + S <-> 


In the limit considered the baryon density is not influenced by the absolute scale of mo but only by 
the ratio /z/ TO o with the mass of the Bino m ^ being the only relevant scale. In particular, in order 
to achieve correct value of VL AB h 2 ~ 0.02 [15] a value g/mo 1 appears favoured. We also notice 

that the factor ^ ) f 2 (''j suggests a suppression of the baryon abundance both for mg <C m^ 


and mg ~ m^. We have for it a maximal value ~ 0.16 for mg/m^ ~ 0.3. We remind however, 
that eq. (33), relies on assumptions valid only in a limited range of the parameter space. We will 
thus postpone a quantitative determination of as function of the MSSM parameters, to the 

next section, once the detailed numerical treatment is considered. Note that the expressions above 
are valid only in the limit in which it is possible to neglect the impact of wash-out processes and 
entropy dilution. 

Wash-out processes guarantee that no baryon asymmetry is created in thermal equilibrium and, 
if they are efficient up to rather late times, they can deplete partially, or even completely, the 
asymmetry created by the decay and the annihilations of the Bino. The main wash-out processes 
are inverse decays of three quarks into a Bino or a Gluino, as well as 2 —> 2 scatterings of the type 


simple expression. This is not fully motivated given the possibility, as shown below, of relativistic or semirelativistic 
decoupling of the Bino. A generalization of the expression is straightforwardly obtained by inserting the cross- 
sections in suitable integrals. 
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uB(G ) didj, diB(G) O udj (and their CP conjugates). In addition one should also consider 
3 —» 3 scatterings of the type udd —> udd , mediated by two scalars and a off-shell gaugino, and, 
similarly, 2 —» 4 scatterings [7]. However these last two kinds of processes have very suppressed 
rates, as m^ 8 , and thus have been neglected in our analysis. A quantitative computation of the 
abundance of baryons including the effects of wash-out processes requires the solution of Boltzmann 
equations and will be discussed in detail in the next section. We can nonetheless distinguish two 
simple limit cases. As mentioned before the baryon asymmetry is mostly generated by the out-of¬ 
equilibrium decay of the Bino with a typical time scale determined by: 


r/j, t ot(*d) « H(x d ) (34) 

By an analogous rule of thumb we can define the scale x w , a , at which wash-out processes become 
inefficient. If -C x w , 0 ,, the baryon asymmetry is produced when the wash-out processes are very 
efficient and, as consequence, it is partially or totally depleted. In the opposite case the baryon 
production occurs, instead, when wash-out processes are not important anymore and, hence, the 
Bino abundance, weighted by the branching ratio of the S-violating processes, is totally converted 
in the baryon abundance. For kinematical reasons, as well as the presence of the strong coupling, 
the most important wash-out processes are the ones related to the gluino, with corresponding rates: 


where z= vn ^-, 

m B 


riD 



Ulf. r, 

~^x 2 K 2 (zx) 

m o 


(35) 


r s 


16a s 
97r 2 


TOn X 


I\i(zx) 
K 2 (zx ) 


K 2 (zx) 


(36) 


describing, respectively, inverse decays udd{udd) —*■ G and 2 —> 2 scatterings, like e.g. ud —> dG. 
These two rates, normalized with H, have been compared with the decay rate of the Bino in fig. Q. 
Here we have considered the following assignments of the parameters: X = 0.1, m^ = 2 TeV, z = 0.5 
and too = 10 55 GeV (left plot) and mo = 10 6 GeV (right plot). In both cases 5 has been kept fixed 
at 10 s GeV. 

For the lowest value of mo the decay of the Bino occurs before the wash-out processes become 
ineffective and we thus expect that at least part of the generated baryon asymmetry is erased. As 
too increases the rates of wash-out processes become more suppressed; the decay rate of the Bino 
is analogously suppressed such that its decay occurs at later times. As shown by the last panel of 
[4] wash-out processes become negligible for mo = 10 6 0 GeV. 

We have therefore a further indication that the efficient production of the baryon abundance requires 
high values of mo, at least > 10 6 GeV. On the other hand we remind that the CP asymmetry e is 
suppressed by TOq 2 and as a consequence a too high mo would lead again to an insufficient amount 
of baryon asymmetry. We thus expect that the correct amount of the baryon asymmetry is achieved 
for a rather definite range of values of mo- 

In addition to the wash-out processes, the produced baryon asymmetry can be as well reduced by 
entropy injection effects. Indeed, as already noticed in 00 , an high enough abundance of the Bino 
can dominate the energy density of the Universe such that its decay is accompanied by a sizable 
entropy injection. We can thus define a dilution factor [5j: 


Ul.85 


l/4^tef.o>B 

*,s /- 


4 = MAX 


r§, tot-^P 1 


(37) 
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FIG. 4: Ratio of the total decay rate of the Bino (red solid curves) and of the two wash-out processes, namely 
inverse decays (blu dashed curves) and 2 —» 2 processes (black dot-dashed curves), related to the gluino, over the 
Hubble expansion rate for the two values of mo reported in the plot and for rrig = 2TeV. In the left plot the de¬ 
cay timescale of the Bino, namely r ~ H, is much lower than the one at which washout-processes become inactive. 
As a consequence the baryon asymmetry is expected to be at least partially erased. For the higher value of mo 
the rates of the wash-out processes are instead below H and the generation of the baryon asymmetry is maximally 
efficient. 


From the discussion above it is evident that the correct generation of the baryon asymmetry depends 
from two absolute scales, being the mass of the Bino m^ and the scalar mass scale mo- The other 
two scales, namely fj, (entering only into the pair annihilation processes into two Higgs) and m.Q 
can be instead determined, as function of, respectively, mo and m^, by requiring the generation of 
the maximal amount of asymmetry. 

The impact on the parameters space of the effects, namely wash-out and entropy dilution, described 
above, as well as the range of validity of the analytical expressions are, qualitatively, described in 
fig. (15]). Here we show the bidimensional plane (m^,mo) for two assignments of the parameter /r, 
namely fi = 10, lOOrno, with tuq = 0.35 rrig and A = 0.3, while the mass of the Wino has been set 
to a much higher scale with respect to the other gauginos in order to decouple possible effects. As 
already argued wash-out processes are active at the lower values of m o- In order to avoid these 
effects we need to require mo to be at least 2-3 orders of magnitude above the scale of the gauginos 
involved in the generation of the baryon asymmetry. The region of impact of wash-out processes 
(green region) substantially corresponds to the scenario in which the typical decay time of the 
Bino, set, by rule of thumb, by the condition r to t = H , is close to the one of freeze-out (yellow 
region). The production of the baryon asymmetry is instead very efficient for much later decay 
times. Entropy dilution effects (light-blue region) occur instead for very high values of mo (and, in 
turn, /i) for which the decoupling of the Bino is relativistic (blue region) while result negligible for 
a production of the baryon asymmetry from out-of-equilibrium decay in the non-relativistic regime. 
We have finally inserted, for reference, the isolines Xd = x^w = ^he re §i° ns the right of 

the curves correspond to a production of the baryon asymmetry before the EW phase transition, 
with its consequent reduction due to sphaleron processes. 

From the discussion above it is thus evident that an optimal production of the baryon asymmetry 
corresponds to a rather definite range of values of mo, m o ~ 10 6 "^' GeV. 
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jL/=10 m Q /7=10 2 m 0 



FIG. 5: Summary plots describing the regions of validity of the analytical estimates in the plane (nig, mo), for 
two values of the ratio /i/mo, namely 10 (left panel) and 100 (right panel). The mass of the Gluino has been set to 
m^ = 0.35 m^ while A = 0.3. The Wino, finally, has been assumed to be very heavy and decoupled from the the¬ 
ory. The green region corresponds to typical decay times of the Bino smaller than the ones at which wash-out pro¬ 
cesses become ineffective. This region substantially overlaps with the yellow region corresponding to decay time, 

Xd, lower than the freeze-out time. In the blue region the Bino decouples while relativistic while the light blue re¬ 
gion indicates sizable amounts of entropy injected at its decay. The production the baryon asymmetry is mostly 
efficient in the white strip outside the regions described above. The dot-dashed magenta line curves correspond to 
the case Xd = x ew (see text for details). In the region below this curves the baryon asymmetry is produced before 
the EW phase transition and it thus depleted by a factor 


B. Production mechanisms for the gravitino DM 


In a supersymmetric scenario there are in general three production mechanisms for the gravitino. 
There is first of all the contribution from thermal scatterings occurring at high temperatures in 
the Early Universe giving a contribution to the relic density sensitive to the gravitino and gaugino 
masses as well as to the reheating temperature after the inflationary phase [50H52] . The contribution 
to the DM relic density is given by ED IM|: 


n ™ h2 ~ (S) (iffiSlev) (1 + S r ) ^1 + l0 S 


(38) 


where y r , k r and 5 r are numerical coefficients defined in ED- In addition we have a contribution 
from the freeze-in mechanism originated by the decays of the superpartners while they are still in 
thermal equilibrium |34| . The expression of the relic density can be written as: 


n 


FIMPi 2 
DM a 


1.09 X 10 27 

3/2 m 3/T ft 

9* 


(39) 


where T, is the decay rate of the i-th superpartner, which can be a gaugino or a scalar, and it 
is given by eq. (221 by substituting the suitable mass, while represent the internal degrees of 
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freedom to the i-th state. Since the decay rate depends on the fifth power of the mass of the 
decaying particle the DM relic density is mainly determined by the decays of the heaviest particles. 
It can be easily seen that it largely exceeds the experimental value because of the high scale of the 
scalars in this setup. The only way out is to impose the condition Xr < Too, in such a way we have 
no equilibrium population of the heaviest states in the early Universe. From now on we will thus 
assume the condition nig < Xr < mo, i.e. a reheating temperature below the mass of the scalars, 
in order to avoid the existence of a thermal population of these particles, but still sensitively above 
the mass scale of the Bino, in order of not affect the generation of its abundance. This requirement 
is not problematic in our scenario since, as shown in the previous subsection, an efficient generation 
of the baryon asymmetry requires a 2-3 orders of magnitude separation between the scales nig and 
Too- The freeze-in relic density thus reduces just to the contribution of the three gauginos which 
can be written as: 



We also remark that our requirement on the reheating temperature implies, as by-product, a sup¬ 
pression of the contribution from thermal scatterings, eq. (38). The DM relic density is then totally 
accounted by the decays of the gauginos. 

We have finally, for the contribution from the out-of-equilibrium decay of the Bino, the SuperWIMP 
contribution: 


n sw 

“dm 


= £, 


m 3 / 2 


Wli 


Br(B —>■ ip 3 /2 + X)Q 


T—> OO 

B 


(41) 


Contrary to the baryon density the only suppression term present is £ s , which accounts for possible 
entropy dilution effects. A similar contribution to eq. (41) originates also from the decays of the 


Wino and the Gluino after they undergone chemical freeze-out. However these two particles have 
a much lower relic density, compared to the Bino, in virtue of their efficient annihilation processes 
and thus the corresponding contribution is negligible. 


Assuming for OjW°° the expression given in eq. (31), we can write: 
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DM 


2.34 x 10 


~ 3 (- 




Y- 
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0’( 


\~Y ^3/2 


Xf.< 


\ 10 3 / 2 ?no/ V 10 6 GeV/ UTeV/ VIGeV/ A{x ( . D .) 


1 + —) (42) 

7ra< 


By comparing eq. (42) and (1401) we notice that the SuperWIMP contribution tends to dominate at 


higher values of mo and of the ratio /i/roo while the freeze-in one becomes more important once 
increasing the masses of the gauginos. In particular the DM relic density can result dominated by 
a heavy Wino, as required, as will be clarified in the next section, by avoiding coannihilation effects 
reducing the baryon abundance. 

If the out-of-equilibrium decay of the Bino is the main source of the DM abundance, the ratio 
Dab/Ddm assumes the simple form, as function of the supersymmetric parameters: 


Da b _ m p ecp 

Ddm m 3/ 2 Br (b -> ^ + X) 



(43) 
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where, in the last line, we have taken the maximal value for jg y r %/ m %J ~ 0.16. Note that 

this ratio is independent of the abundance of the decaying Bino and that, interestingly, the correct 
ratio between the two relic densities is achieved, for a Bino at the TeV scale, when the gravitino 
mass is of the same order as the mass of the proton. Unfortunately, as clarified by the numerical 
treatment in the next section, the requirement of the correct abundance of the Bino, mandatory for 
the matching of the individual quantities with their observed values, will point towards sensitively 
higher masses for the Bino and the gravitino. 


V. NUMERICAL ANALYSIS 


A. Boltzmann equations 


The generation of the baryon asymmetry and of the DM, including additional effects like wash-out 
and coannihilations, in the scenario under study, can be traced through a system of five coupled 
Boltzmann equations. The first three describe the evolution of the yields, namely Y = n/s, of the 
three gauginos: 
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(46) 


In each equation the first row represents B-violating decay and single annihilation processes. The 
second to the fourth lines represent coannihilation and pair annihilation processes. The remaining 
lines, apart the last, give rise to transition processes, either decays or scatterings, between gauginos. 
The last line in each equation represents finally the production of the gravitino. These last decay 
terms are proportional only to the yields of the gauginos since we assume that the initial gravitino 
abundance is negligible and remains low enough to neglect inverse decay processes. Under this 
assumption the equation for the gravitino abundance assumes a rather simple form: 

7 i = 7Y r ( i ^ + x ) Y * * = < 47 > 

x 


For simplicity we are neglecting the possibility that the Bino dominates the energy density of the 
Universe since, as already argued in the previous section and further confirmed by the results 
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presented below, this occurs in a region of the parameter space of marginal relevance. As a conse¬ 
quence, the expression of the Hubble expansion parameter is the one typical of radiation domination, 

m 2 

H ss I. 665 * m^ x ~ 2 - or( ler to properly account entropy injection effects it should be modified 
similarly to what proposed e.g in [35H38] . We have finally the equation for the baryon asymmetry 
which is casted as an equation for Yab-l in order to get rid of the effects of the sphalerons: 
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The first row represents the source terms associated to the B-violating decays of the Bino and, 
as already mentioned, to the scatterings of both Binos and Gluinos. CPT invariance imposes a 
relation between the asymmetries generated by Binos and Gluinos mmm-- 


(Aav) s Y~ q = - (A av) 6 Y^ 


(49) 


In general we could expect analogous source terms associated to decay and scattering processes with 
Wino initial state. As already discussed in the previous section and shown in an explicit example 
below, the Wino is always kept very close to thermal equilibrium by its efficient interactions and thus 
contribute to a negligible amount to the generation of the baryon asymmetry. The last two rows 
describe instead the wash-out processes related to inverse decays and to the CP even component 
of the baryon number violating 2 —f 2 scattering of both Binos and Gluinos. The equation for 
the baryon asymmetry depends as well on the chemical potentials fif= u ,d,s,c,b,t of the right-handed 
quarks. These chemical potential can be expressed in terms to of the B — L abundance [3?J . We 
have in reality different relations between the chemical potentials and B — L according on whether 
the temperature lies above or below the one of the EW phase transition. Since, as will be discussed 
below, the production of the baryon asymmetry can occur, according the values of the relevant 
parameters, both above and below this critical temperature, we have employed, similarly to what 
done in [5], a two step solution of the system eq.(44)-(|48|. We have first solved the system with 


initial conditions Y^(x « 1) = Y? q (x -C 1), F v p( x <C lj = Y^(x -C 1), Yq(x « 1) = Y? q (x < 1 ) 
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and Yab-l(x -C 1) = 0 and: 


Hu — He — Ht — 


10 Yab-lS 
79 T 2 


5 Hd — Hs — Hb — 


38 Yab-ls 
79 T 2 


(50) 


unt il x = nig/T ew- Below the EW phase transition the sphalerons freeze-out and we can replace 
eq.(48) with an equation for just Yab with initial condition, set at Tew ; Y/\b = MYb-l and: 


Hu — \ L + B 
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(51) 


where: 


12tt 2 £*,s 

-d = -—- rrigXYAB 


L = 


45 
127r 2 ff 4 ,g 
45 


rrigxY L (TEw) 


(52) 


The structure of the system makes evident the tight relation, already envisaged in the analytical 
treatment, between the generation of the CP asymmetry, the abundance of the Bino and the wash¬ 
out processes. The baryon asymmetry is originated by the decays (and annihilations) of the Bino 
into SM fermions. Inverse decay and scatterings, as well as analogous processes involving the Gluino 
(the processes related to the Gluino are generated by the same diagrams and then the rate differ 
only by the couplings and by the Gluino distribution) are responsible of the wash-out of the baryon 
asymmetry. Single RPV annihilations of the Bino can, finally, be the dominant contribution in 
determining its abundance. Any variation in the CP asymmetry, as originated, e.g. by flavor or 
left-right mixing effects, is reflected also in this last rates. The optimal production of the baryon 
asymmetry is thus achieved once a balance is found between a large enough CP-asymmetry and 
not excessive depletion of the Bino abundance, or excessive efficient wash-out. 

The system has been solved for several assignments of the relevant parameters. Differently from 
the analytical treatment we have considered also the cases in which the relevant processes are 
mediated by left/right-handed top squarks, as well as generic effects of flavor violation in the right- 
handed down squark sector by assigning arbitrary entries and CP violating phases to the matrix 
Tg and taking non-degenerate squark masses. In both these two cases we have found no sensitive 
variations with respect to the flavor universal scenario. In the case of left-right mixing in the top 
squark sector this is due to the fact that the (very moderate) enhancement of the CP asymmetry 
is actually compensated by the presence of the electric charge Q u of the up-quarks in the couplings 
of the Bino which translates into an overall increase by a factor 4 of the annihilation rates of the 
Bino (it can be easily seen from the analytical expressions which, on the contrary, the value of 
the CP-asymmetry is insensitive to this quantity.). In the flavor violating case instead the small 
variation in the total CP-asymmetry is due to the already mentioned GIM suppression. For this 
reason we will discuss our results in the same flavor universal limit of the analytical treatment in 
order to profit of the more limited set of parameters, being (A, nig, tUq, mo, /r). In all cases 
it has been found that the dominant contribution to the baryon asymmetry is originated by the 
out-of-equilibrium decays of the Bino. Several examples of numerical solutions will be illustrated in 
the next subsections. We will first of all show quantitatively the effects of wash-out processes and 
the impact of the Wino in the generation of the baryon asymmetry. We will then determine the 
regions of the parameter space which provide the correct baryon abundance and the correct DM 
relic density. 
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B. Effects of coannihilations and wash-out 


We show in the following some examples of numerical solutions of the system of Boltzmann equations 
highlighting in particular the impact of coannihilations and wash-out effects. 



FIG. 6: Evolution of the abundance of the yield Yg (solid lines) and of the baryon abundance (dashed lines) with 
x = m,g/T for m,g = 2TeV,m ( j = 1 TeV, A = 0.1, // = 10 8 GeV and four values of mo ranging from 10 5 ° to 
10 7 GeV reported in the plot. 


Fig. © shows the Bino (solid lines) and the baryon yield (dashed lines) for several values 
of too, ranging from 10 5 5 to 10' GeV, and with the following assignment for the remaining 
parameters:???^ = 2 TeV,= ITeV, A = 0.1, /i = 10 8 GeV. For the lowest values of mo we have 
a low baryon abundance as consequence of the suppressed abundance of the Bino, whose yield 
remains close to the equilibrium distribution until late times. Moreover the baryon abundance is 
almost completely depleted for mo = 10 5 GeV since for this values of the scalar mass scale the Bino 
decays before that wash-out processes become ineffective. The baryon abundance is maximal in 
the intermediate mass range, order of 10 6 ' 5 GeV, where the Bino features a rather early decoupling 
and it is long-lived enough to evade the wash-out regime. The baryon density then decreases again 
at higher masses when the Bino gets close the relativistic decoupling. Indeed its relic abundance 
is poorly sensitive to the increase of mo while the CP asymmetry ecp still features a sensitive 
suppression. This result justifies our choice to neglect eventual deviations from standard cosmology 
in the numerical system. Indeed entropy production occurs in the very high mo region which is not 
relevant for our analysis since we expect a suppressed asymmetry. 

Fig m shows the evolution of the abundances of the Bino and of the baryon density, compared 
with the ones of the other two gauginos. As evident these two species tend to remain in thermal 
equilibrium (up to their decay) during the whole phase of generation of the baryon asymmetry. 
The four panels of fig. (j7| differ in the assignments of the mass of the Wino, considered to be 
both below and above the mass of the Bino. The Wino has a profound impact in the generation 
of the baryon asymmetry. The case of a light Wino is, in particular, disfavored. Indeed in such a 
case, coannihilation effects turn to be very strong, keeping the Bino very close to the equilibrium 
distribution up to late times, with consequent suppression of the baryon abundance. Contrary to 
conventional WIMP coannihilation scenarios, the Bino abundance results altered even for sizable 
mass splitting with the Wino. This is consequence of the strong suppression of the Bino annihilation 
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FIG. 7: Evolution of the Yields of the Bino (red solid line), the Gluino (black solid line), Wino (blue solid line) 
and of the baryons (green solid line), for three assignments of the mass of the Wino. For reference has been re¬ 
ported as well the equilibrium distribution of the Bino (red dashed line). 


rates. This last effect is better evidenced in fig. (|8| where even higher values of the ratio have 
been considered. In order to maximize the production of the baryon asymmetry we need to invoke 
a strong hierarchy between the mass of the Bino and the one of the Wino, at least > 5. As 


argued in the previous subsection, such a heavy Wino can copiously produce DM through freeze-in. 
This issue can be possibly avoided by requiring a very heavy Wino with > Tr or by checking 
that it is light enough to avoid overclosure as given by imposing the condition VL^ M h 2 < 0.1, i.e. 
from eq. (401 


Tow < 362 TeV 


/ W'3/2 

Vl TeV 


1/3 


(53) 


In the next subsection we will focus on the case in which the decays of the Bino are the primary 
source of DM production and we will thus assume, for simplicity that the mass of the Wino is above 
the reheating temperature. 
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FIG. 8: Evolution of the abundance of Yg (solid lines) and of the baryon abundance (dashed lines) for a definite 
assignment of (A, nig, mg, mo, /i and varying m^/m,g), as reported in the plot. 


C. Results 


We will illustrate below the regions of the parameter space accounting for the experimentally favored 
values for the baryon and DM abundances. In our setup the baryon asymmetry depends on five 
parameters, namely the mass of the Bino mg, the mass of the gluino mg, the heavy scales mo and 
/z, and the RPV coupling A. The DM relic density depends on two additional parameters, the mass 
of the gravitino 7713/2 and, possibly, the mass of the remaining gaugino m Regarding this latter 
parameter, as already discussed, a value close to the masses of the other gauginos is disfavoured 
by the correct baryon asymmetry. We will, from now on, implicitly assume, for simplicity, that the 
mass of the Wino is decoupled from the relevant phenomenology, i.e. m^ v > Tr. 

As discussed above, the baryon density is the most difficult quantity to accommodate. We will thus 
determine it in the bidimensional plane (nig, mo) after having identified an optimal assignation for 
the remaining parameters. The correct DM abundance can be determined accordingly by a suitable 
choice of the mass of the gravitino. 

Fig. 0 shows the evolution of the yields of the baryons and of the DM as the parameters mg 
(left panel) and A (right panel) are varied, while keeping fixed the others. As evident, in the case 
of both quantities there is a non-trivial interplay in the determination of the baryon abundance. 
A mass of the gluino very close to the one of the Bino determines a huge suppression of ecp (in 

eq. 24 /1 -C 1) while in the opposite scenario, i.e. mg/?77/j <C 1 the baryon abundance is 

analogously suppressed by the factor mg/m/j in ecp and, more important, the wash-out processes 
are efficient up to very late time scales, with respect to the one of decay of the Bino, substantially 
depleting the created asymmetry. This is then maximal for r ^- L ~ 0.3 — 0.6. For such values there 
is still a sizable kinematic suppression of the B -violating decay of the Bino as well as its abundance 
due to the effect of the coannihilations with the gluino as well as the single annihilations into a 
Gluino final state. A similar situation occurs also for the A coupling, with a suppression of the 
baryon abundance both for A ~ 1 and for A < 1. The behaviour at high values of A is motivated 
by the fact that ecp is independent of such coupling in this regime (see eq. 24 ). As a consequence 


the main effect is the increase of the rate of the single Bino annihilations, influencing both the 
Bino abundance and, directly, also the one of the baryons, through an enhancement of wash-out 
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FIG. 9: Left panel: Bino (solid lines) and baryon yields (dashed lines) for a fixed assignment of mo and /x, 
reported on the plot, for A = 0.3 and for four values of itlq ranging from, 0.05m^ and 0.9m g. Right panel: The 
same as left panel but with rriQ = 0.35and A varying between 0.05 and 0.9. 


effects. In the regime of very low A the dominant effect is the suppression of the branching fraction 
of ^-violating decays since the abundance of the Bino is controlled by the annihilations involving 
the Gluino as well as the pair annihilation processes. The optimal range for the A parameter is, 
again, the intermediate range A ~ 0.3 — 0.6. 



m ~[GeV] 


m g[GeV] 


FIG. 10: Contours of values of the baryon abundance Yab in the plane (m^,mo). The /x parameter has been 
set to 100mo- In the left panel the Bino mass have been varied in the range 1 — 10 TeV while in the right 
panel higher masses of the Bino, namely 20-100 TeV, have been considered. For both plots we have considered 
m^/m^ = 0.4 and A = 0.4. 
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Fig. (101 reports the isocontours of the baryon abundance Ya/j in the bidimensional plane (m^j, too) 
with m.Q and A fixed, according the discussion above, to, respectively, 0.4 rrig and 0.4. As already 
argued in our analytical study the correct order of magnitude is achieved only for a rather restricted 
range of values of too- Above this region there is an excessive suppression of ecp while below the 
baryon abundance is erased by wash-out processes. The correct value of the baryon abundance is 
achieved for a rather heavy Bino, with mass ~ 70TeV, and Too ~ 10 7 ° GeV (by varying the ratio 
TOq/to g and A within the range indicated above it is possible to lower to approximately 50 TeV the 
minimal viable Bino mass). This is due to the suppression, direct and indirect, of the Bino density 
coming from the presence of a rather close-in-mass gluino, which requires high values of the scales 
rrig and too to be compensated. As can be seen from the second panel of fig. (101, the correct 
baryon abundance can be achieved also for > 100 TeV. However the consequent increase of the 
scale mo would create tension with the determination of the Higgs mass gnmmm- 
As already mentioned we have assumed throughout this work the tan/3 —> 1. As shown in m it 
is possible to have the correct value of the Higgs mass also for tan/3 > 50 and too = |uL We have 
solved the Boltzmann system and found analogous contours as the ones shown in fig. ( fl0| ) also for 
/< = too- In this case BB —► HH* annihilations play no relevant role and we can thus reduce the 
number of free parameters, although the general results remain substantially unchanged. 


jL/=1 00/77 q 



m ~[GeV] 


FIG. 11: Isocontours of the baryon and DM Yields. The red band represents the value Yab = (0.86 ± 0.01) X 
10“ 10 determined by CMB measurments ESI. The blue dashed lines represent the extrema 0.4 x 10 10 and 
0.9 X 10 —10 determined by BBN 141| . The black, gray and purple dot-dashed lines represent the isocontours of 
the correct DM relic density for the reported values of the mass of the gravitino. 


The baryon abundance is finally compared with the one of DM in fig. (11). Here we have reported 
the experimentally favored value Tab = (0.86 ± 0.01) x 10 -11 and confronted it with isocontours 
of the correct DM relic density for some values of the gravitino mass. As we see the correct match 
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between the two abundances occurs for a mass of the gravitino between, approximatively, 100 GeV 
and 3 TeV. A lower mass of the gravitino is achieved if a wider range of variation, like the one 
shown in the figure based on BBN measurements ED, is allowed. 


m 3B =100,m~ =7.2*10 4 , m 0 =10 7 5 , /J=10 95 


m 3 ,2=3000,m~=9.7*10 4 , m 0 =10 78 , /J=10 98 



FIG. 12: Two benchmarks featuring the correct ratio between the DM and baryon abundances, as well as the 
correct agreement of the individual quantities with the experimental determination. In both cases the DM and the 
baryon asymmetry are produced by the out-of-equilibrium decay of a semi-relativistically decoupling Bino. 


As shown in fig. the optimal benchmarks highlighted in fig. ( |11| ) (namely the gray and purple 
lines), correspond to a contemporary production of the DM and of the baryon asymmetry from 
the out-of-equilibrium decay of the Bino with the latter featuring a substantially semi-relativistic 
decoupling. We also notice that the yield Y 3 / 2 of the DM is sensitively lower than the one of the 
baryons but it is compensated by the much higher mass, with respect to the one of the proton, such 
that the relic density results bigger, as expected. 

The result obtained is sensitively different with respect to the scenario proposed in na, consisting in 
accommodating the correct value of Dab/Udm through similar values of ecp and of the branching 
fraction of the mother particle into DM and, accordingly, similar values of the gravitino and proton 
masses. The reason of this resides in the more accurate determination of the abundance of the 
Bino as well as of ecp, in particular the inclusion of the kinematic functions /i, 2 , and especially 
in the impact of wash-out processes. The combination of these effects leads to the need of a 
heavier supersymmetric spectrum to achieve successful baryogenesis and therefore a substantially 
heavier gravitino is required to compensate the suppression of the Bino branching fraction into DM. 


Indeed, the quantities ecp and BR yB —> ip 3/2 + A'J differ in the allowed window by a few orders of 
magnitude and the similarity between DM and baryon densities cannot be directly related to their 
near-equality. At the same time we remark that we could achieve a viable scenario in which the 
correct amounts of baryon asymmetry and of DM are contemporary produced by the decay of the 
Bino at relatively low value of the reheating temperature. Our predictions depend, apart a single 
assumption on the cosmological history of the Universe, on the masses of the superpartners, and we 
were able to identify rather definite ranges for the supersymmetric particles masses, in particular 
the gauginos. 

Before concluding this section we just mention that an extensive exploration of the space of the 
parameters (in particular the flavor-mixing matrices T^) involved in the generation of the baryon 
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and DM densities is rather complex. Although we have found, in our analysis, that, as discussed 
above, there is no particular loss of generality in assuming contribution from only degenerate d- 
squarks to the relevant processes, we cannot completely exclude the presence of configurations, 
possibly involving also contributions from u-squarks, leading to 0 ( 1 ) variations, with respect to the 
results presented here, might occur. We notice in particular from fig. (10 ) that an enhancement of 
a factor 2 — 3 of the CP asymmetry, not compensated by annihilation or wash-out effects, would 
allow for smaller masses of the Bino and the Gluino, below 10 TeV, with the latter possible even 
lying in the LHC production range. 


VI. DETECTION PROSPECTS 


In this section we will briefly investigate possible experimental signatures of this scenario and bounds 
associated to them. The main experimental signature of our scenario is the Indirect Detection of 
the decays of the gravitino DM. Indeed due to the RPV coupling A , the gravitino has a three-body 
decay into SM quarks, possibly leading to signatures in the antiproton spectrum, with a rate m-- 
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(54) 


with N c being the number of channels giving a lifetime: 
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Interestingly, for values of m 3/2 and mo of, respectively, ITeV and 10 7 ' 5 GeV, which provide the 
correct fit of the DM and baryon abundances, a DM lifetime of approximately 10 28 s is achieved, 
which is exactly of the order of the current AMS-02 sensitivity in the antiproton channel ii 
and thus allows to test in the very next future our scenario. The decay into quarks of the gravitino 
can give rise as well to a sizable signal in 7 -rays. Similar sensitivities, to the one discussed for 
AMS-02, are expected for 7 -ray detectors like H.E.S.S. and CTA [45]. The heavy supersymmetric 
spectrum does not offer, instead, very promising prospects for collider detection. The scenario 
proposed requires possibly a supersymmetric spectrum beyond the kinematical reach of LHC while 
the Gluino NLSP could be within the expected reach of a 100 TeV collider |46tf4%] . As mentioned 
above, within the factor one uncertainty of our computations, we cannot exclude the possibility of 
having a slight enhancement of the CP violating parameter ecp allowing for viable baryogenesis 
and DM production in regions of parameter space with a lighter supersymmetric spectrum. In case 
of a mass of < 10 TeV it would be possible to observe the Gluino NLSP at the LHC. In our 
scenario its main decay processes would be mediated by the RPV coupling A with typical decay 
length: 


2.75 / A \ 2 / too \ 4 / \ -5 

CT ~ 9 ~ ~N2 m \(u) V 10 7 GeV/ \2TW/ 


(56) 


corresponding to displaced vertices or, most probably, a detector stable state. The prospects of 
detection can be inferred using the techniques discussed, for example, in [49l l5T| . The detection of 
EW gauginos requires anyway next future, higher center of mass energy facilities | I6I11<8| . 

We have not discussed here any particular flavour structure of the RPV couplings A, but in principle 
those couplings can also contribute to flavour-violating neutral current processes, as well as B- 
violating processes (other then proton decay) like neutron-antineutron oscillations and A B = 2 
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transitions (see e.g. [TB] for an extensive discussion.). Since in our case we need a very large scale 
for the scalar quark partners, even if we need a large coupling, those rates remain well below the 
present limits and will be difficult to reach also in the future. 


VII. CONCLUSIONS 

We have presented a systematic approach for determining the contemporary production of the DM 
and the baryon asymmetry from the out-of-equilibrium decay of the same mother particle in a 
MSSM framework. These two quantities have been numerically determined through the solution 
of a system of coupled Boltzmann equations, accurately computing the abundance of the decaying 
state and taking, in particular, in account of the impact of wash-out processes. We have supported, 
whenever possible, our numerical results with analytical estimates. We have determined the ranges 
of the values of the relevant supersymmetric parameters which allow for an efficient production 
of the DM and baryon abundances. In the most simple (but rather general) limit of only right- 
handed d-squarks involved in the generation of the baryon and DM densities, the observed ranges 
for those quantities are met for a value of the mass of the decaying Bino of 50-100 TeV, a mass 
of the gluino NLSP of 20-50 TeV and a mass for the DM gravitino LSP between 100 GeV and a 
few TeV, and all the other supersymmetric particles above the scale of 10 7 GeV and not present 
in the primordial Universe because of the assumption of a lower reheating temperature, in order 
to avoid the overproduction of the DM. But note that a slight increase of the ecp parameter by a 
factor 2 or so, due to the presence e.g. of intermediate up-squarks, if not completely compensate 
by an increase in the wash-out processes or the Bino annihilation rate, could allow to reduce the 
supersymmetric masses by a factor of a few. 

We find moreover that the similarity of the DM and baryon densities cannot be explained by the 
relation ecp ~ BR (j3 —> DM + anything^ and the gravitino mass has to be tuned to give the 

correct DM abundance. Nevertheless the common generation of baryon and DM density from the 
Bino neutralino after freeze-out can work and provide the right abundances for large values of the 
RPV coupling and in cosmologies with low reheating temperature. 

The very heavy supersymmetric spectrum does not offer promising detection prospects at the LHC, 
but the Gluino LSP could be within the reach of a 100 TeV collider. On the other hand a very 
promising signal in the near future might come from the decay of the gravitino whose lifetime can 
be within the present sensitivity of AMS-02 and gamma-ray detectors. 
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Appendix A: General expressions for the CP asymmetries 

In this appendix we will provide some general expressions, including the dependence on the flavor 
matrices T^’^ of the CP-asymmetry possibly originating the baryon abundance in our setup. As 
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shown by our analysis the baryon asymmetry is mainly generated by the Bino decay. For sim¬ 
plicity we will just focus on the CP asymmetry in the decay processes. The computation can be 
actually straightforwardly extended to the 2 —> 2 scatterings since the corresponding rates are 
cross-symmetric to the decay ones. 

As already mentioned in the main text the CP asymmetry is defined as: 


ecp 


( B Ukdidk^j - P (^B -> Ukdidkj 
{jd Ukdidkj + r ^ B — > Ukdidkj 


(Al) 


The total CP asymmetry is given by the sum of the single CP-asymmetries in the different channels, 
weighted by the correspondent branching rations. As well known a non-zero CP asymmetry requires 
the interference of tree and loop level contributions. In the scenario under considerations the relevant 
decays of the Bino are three-body processes in three SM quarks (un up-type quark and two d-type 
quarks) or a gluino and a SM quark pair (the two body decay into DM is irrelevant for the generation 
of the baryon asymmetry). These processes are mediated by down and up type squarks. 

For the case of the decay into only SM states we have that: 
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In the case that the squark matrices are diagonal and the down-type quark are sensitively lighter 
than up-type quarks we have that: 

r (b ->• u k did k + B ->• Ukdidk) = |Afc £f^ (A 3 ) 

V / 12871^ rriQ 

The other relevant rate at the tree level is the one of decay into Gluino. It is given by: 
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where: 
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We can now move to compute the CP asymmetry. Parametrizing the loop amplitude as Ai 00 p .Fi 00 p 
where A is a numerical coefficient depending on the coupling and the effective CP-phase while F is 
a suitable loop integral we have that: 

AP = T [b -» u k did]}j - T (j3 -)• Uhdidkj = 4:Im(Al ree A loop )Im(F\ oop ) (A 6 ) 

In order to properly identify the different contributions to the CP-asymmetry, in particular in 
relation to the flavor structure, we have performed our computation using the fully Supersymmetric 
Lagrangian, rather than dbl, and performed at the end of the computations the limit rriq m^, mg. 
In the most general case the CP asymmetry originates from a combination of several tree level and 
loop diagrams. We have first of all the diagrams in which only d-type squarks are exchanged. These 
have been reported in fig. (|T|) and[ 2 ]and consist in two tree-level diagrams, labeled as T 1 and T 2 and 
4 loop diagrams, labeled LI, L2, L3 and L4, in which two d-squarks are exchanged. We have then a 
tree level diagram and loop diagram with the same topology as, respectively, T1 and L3, but with 
up-type squarks exchanged. We have finally diagrams with, again the same topology, as T1-T4, but 
with exchange of one up-squark and one d-squark. As already argued all the possible topologies of 
diagrams are already accounted by the case of exchange of only down-type quarks. For this reason 
we will focus on this case since the remaining contribution can be straightforwardly obtained from 
the expressions presented. We show below the values of the decay asymmetry originating from the 
combinations T1L1, T2L1, T3L1 and T3L2. The other combinations are obtained from this by 
exchanging the flavor indices of the two final state d-quarks. 
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The imaginary part of this integral, as well as the others appearing in the expressions below, can 
be computed with the Cutkosky formalism putting the internal quark (with four-momentum l—pk) 
and gluino (with four-momentum l + Pj) on-shell. The integration over the phase space leads to: 
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2. T2L1 
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where: 
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The integration over the phase space leads to: 
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This and the previous expression depend on the modulus square of the effective gauge couplings g G 
and g^. A non null CP asymmetry from the corresponding diagrams might arise only in presence 
of flavor violation and CP violating phases in the T^ L matrices. 
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3. T1L3 


The contribution to the CP-asymmetry associated to this topology is: 
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where C 2 = 4/3 is a color factor arising from the coupling of the Gluino and Xi = 2 Ei/m^. As 
we notice the effective couplings of the Bino appear with the same conjugation opposite to the 
one of the coupling of the Gluino. The Maiorana phases are thus enough to originate a non null 
CP-asymmetry with effective phase <f> = 2{(f> g — 4 >b)- In this case however a suppression factor 
TOq/to^ is however present. I 4 and 1 $ are the loop integrals defined as: 
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Performing as well the integrations over the phase space it is possible to obtain, in the limit m.q^ 
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4. T2L3 


This contribution is computed analogously to the previous one and results to be: 
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The loop integrals are now: 
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Performing the suitable integrations one obtains: 
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Appendix B: General expressions annihilation cross-sections 


We provide in this appendix the complete expressions of the scattering cross-section of the Bino, 
including the flavor structure as well as the exchange of it-type squarks, involving SM fermions. 
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These are: 
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As evident there is a tight relation between the flavor structure of these expressions and the one 
of the decay rates shown in the previous appendix. In the limit mj <C mj and absence of flavor 
violation and left-right mixing, taking also for simplicity rriQ = 0, we obtain: 
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